Asymmetric exchange between electron spins in coupled semiconductor quantum dots 
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We obtain a microscopic description of the interaction between electron spins in bulk semiconduc- 
tors and in pairs of semiconductor quantum dots. Treating the kp band mixing and the Coulomb 
interaction on the same footing, we obtain in the third order an asymmetric contribution to the ex- 
change interaction arising from the coupling between the spin of one electron and the relative orbital 
motion of the other. This contribution does not depend on the inversion asymmetry of the crystal 
and does not conserve the total spin. We find that this contribution is ~10 -3 of the isotropic ex- 
change, and is of interest in quantum information. Detailed evaluations of the asymmetric exchange 
are given for several quantum dot systems. 



O 

o : 

<N 
> 
O 

in 

O 

in 
o 



a 

o 
o 



Spin in quantum dots (QDs) are attractive candidates 
for qubits in quantum information in part because of 
their long coherence times Q. Controllable coupling be- 
tween these spins is an essential requirement for two qubit 
quantum gates, and it has been the focus of much re- 
cent research 0-0. Typically the spin-spin coupling is 
dominated by the isotropic exchange interaction JSrS2, 
which arises from the Coulomb interaction and the Pauli 
principle. The isotropy of this part of the exchange im- 
plies that the total spin is conserved. This is important 
in gate operations, which involve pulsing J(t) 

In general, spin-orbit coupling in solids also gives rise 
to terms that are asymmetric between the spins, and 
which do not conserve total spin. These terms can 
cause loss of fidelity in gate operations involving two 
spins. To date, the asymmetric terms that have been 
discussed involve spin-orbit coupling of individual elec- 
trons in the conduction band 0, 0- For systems with 
bulk inversion asymmetry, this coupling has the Dres- 
selhaus form @, Ej. I n addition, heterostructure asym- 
metry can introduce the so-called Rashba coupling flOj . 
These two couplings can give an additional interaction 
between the spins of the Dzyaloshinskii-Moryia (DM) 
form /3 so -(S 1 xS 2 )/ v /^ 2 o + ^ 2 E3 where (3 so is lin- 
ear in the spin-orbit coupling strength ja, • ^ n general 
the dephasing caused by the DM contribution cannot be 
totally eliminated, but suggestions have been made for 
gate pulse shapings that eliminate it to first order pjj . 

Here we give a different contribution to the spin- 
spin coupling of two electrons in III-V semiconductors 
that does requires neither bulk inversion asymmetry or 
Rashba coupling. It arises from the Coulomb interaction 
between two electrons and from the conduction-valence 
band mixing. It is similar to the exchange interaction 
between excitons 14]. We describe this asymmetric elec- 
tron spin coupling in bulk materials, and we make de- 
tailed evaluations for spins in several coupled QD sys- 
tems. 

We use a 8 x 8 Kane Hamiltonian to represent the band 
structure of III-V semiconductors Q, with a gap be- 
tween an s-like conduction band and p-like valence bands. 
The band coupling in this Hamiltonian is obtained in k-p 



effective mass perturbation theory, where p is the rel- 
evant momentum operator. The standard parameters 
in this approach are the band gap E g , the energy of 
the split-off band A, and the valence-conduction band 
coupling P=(h/mo)(S\p x \X), where \S) and \X) are the 
Bloch states of the conduction and valence bands |l5| . 

We consider two electrons of relative coordinate 
r=i*i— Y2 in a semiconductor of dielectric permittivity 
ft. The unperturbed two-particle Hamiltonian in the 
conduction band is i^°) = (p^+P2)/2mo. The two band 
mixing terms ki-pi, k2P2, and the Coulomb interaction 
Uc=e 2 /hzr between the electrons are treated on equal 
footing. The leading part of the spin-dependent inter- 
action arises in third order from first order contribu- 
tions of Uc and of each of the band mixing terms ki-pi, 
k2-p2- The spin-dependent part of the Coulomb interac- 
tion arises from the coupling between the electron spins 
and their relative motion: 



(3) 



e 2 2P 2 A(2E g + A) 
~ft~3£f 
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The coefficient of the square bracket in Eq.fQJ is the 
coupling constant, which we will call C. For GaAs 
(7=5.7 meVnm 3 , and for In As C=10 meVnm 3 . The re- 
maining spin-independent terms in this order of the the- 
ory contribute to the isotropic exchange J. These include 
a 'local contact' term of the form 
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Coupling terms similar to Hs° J and are known for 
two electrons in the free space 16] where they arise from 
electron-positron band mixing and Coulomb interaction. 
In the present case the effect is stronger because the en- 
ergy gap Eg between electrons and holes in a crystal is 
much smaller than the energy gap m^c 2 between elec- 
trons and positrons. Also, the symmetry of interaction 
Hg 3 ^ [Eq.fDl is different from that between electrons in 
free space 16], because the valence bands and I\ have 
a symmetry different from the s-symmetry of positrons. 
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We note here that for itinerant electrons the interac- 
tion # s (3) [Eq.flJ] flips a spin in an electron-electron col- 
lision, providing a new mechanism for the relaxation of 
spin polarization in addition to other known mechanisms 
[l7j . We evaluated the spin dephasing time for a 2D 
electron gas in GaAs quantum wells and found that it 
is of the order ^1 ns for He temperature and of the or- 
der ~1 ps for room temperatures. This suggests that the 
interaction can be important in spin transport in 

low-dimensional structures, such as quantum wells [l8|. 

Here we will primarily address systems with spins on 
two centers, such as two QDs or two charged donors, 

(3) 

where Hs gives rise to a DM coupling. The two elec- 
trons can be in a singlet state with total spin 5=0 or in a 

(3) 

triplet state with total spin 5=1. Hs has non-zero ma- 
trix elements between states of different total spin. We 
take the lowest singlet and triplet states to be separated 
by an exchange energy J. Then in the Hilbert space of 
these two states the Hamiltonian can be written in the 
form 

H = JSi-S 2 + i/3-(Si-S 2 ) , (3) 

where i/3=—C/h 2 {So\(r/r 3 ) xpi |T ) is the matrix ele- 
ment between the lowest orbital singlet |5o) an d the 
triplet | To). A block diagonal form equivalent to Eq.0 
can be obtained by an orthogonal transformation in spin 
space, giving the so called 'twisted spin representation' 

m 

H = Jcos0Si-S 2 + 2J^sin0 (n ■ Si)(n ■ S 2 ) + 

+ Jsin^n- (Si x S 2 ) , (4) 

where the spin-0 and spin-1 states are mixed by 
the operator exp [i(0/2)n • (Si — S 2 )], with n=(3//3 and 
</>=arctan(/3/ J). In Eq.Q the first two terms are 
the symmetric isotropic and the symmetric anisotropic 
Heisenberg terms. The last term is the asymmetric ex- 
change in the DM form /3-(Si xS 2 )/ v / /3 2 + J 2 . The con- 
tributions from the Dresselhaus and Rashba couplings to 
the anisotropic and asymmetric exchange were given in 
the form of Eq.QJ in Ref.@. 

We consider the lowest single-particle states \cp±) from 
the two confining centers. In the absence of spin-orbit 
coupling the two-particle wavefunctions can be written 
as products of orbital states and spin states. The axial 
vector (3 is nonzero only for inversion- asymmetric con- 
fining potentials. For such systems we take into account 
the possibility that both electrons occupy the same site, 
therefore we use the Hund-Mulliken description of the 
two-particle orbital states [4]. We obtain the ground 
state by diagonalizing analytically the tunnelling and the 
Coulomb and local contact interactions. We then focus 
on the Hilbert space determined by the lowest singlet 
(5=0) and triplet states (5=1) in which the Hamiltonian 
is given by Eq.JSJ). The importance of the resulting DM 
asymmetric exchange depends on the ratio tan (/>=/?/ J 



between the coefficients of the third and the first terms 
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FIG. 1: (Color online) (a) Dependence of the asymmetric 
exchange /3/J on the separation d along the growth axis be- 
tween two cylindrical vertically-coupled dots with a lateral 
offset b—4 nm, for a smaller dot size a_=5 nm and differences 
between radii of dots Aa=a+-a_£[5, 6] nm. (b) dependence 
of asymmetric exchange on the lateral offset b at a separa- 
tion d—10 nm for a_=5 nm and varying Aa. The insets in 
(a/b) are lateral/ vertical views of the geometry, (c) Asym- 
metric exchange [3/J obtained here (solid line) compared with 
the asymmetric exchange [3d /J from Dresselhaus coupling 
(dashed line) as a function of separation d, for Aa=0.22 nm, 
a_=5 nm, b—A nm. 

Consider first the ratio (3/J for electrons confined 
by two vertically coupled QDs, such as to those in 
InAs/GaAs systems [2(|. These dots generally have dif- 
ferent sizes and can have shape asymmetries as well |2l| . 
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The band offset between GaAs and InAs eV) is 

larger than the quantization energy in the lateral direc- 
tion (hw~20 meV) by more than an order of magnitude. 
This allows us to decouple the vertical and lateral de- 
grees of freedom. We describe these dots by the poten- 
tial offsets along the growth axis z and parabolic poten- 
tials in lateral directions x and y, which in general can 
be anisotropic. These lateral directions are independent 
of the crystal axis. The potentials and wavefunctions 
then can be written straightforwardly 0,Hll- Thus, the 
asymmetry studied here results from differences in the 
lateral sizes and shapes of the QDs. 

The dependence of the ratio (3/ J of the asymmetric 
exchange to the symmetric exchange on the separation d 
between the two cylindrical dots with a fixed lateral off- 
set b is shown in FigEJa). Its dependence on the lateral 
offset b for two dots at a fixed d is shown in FiglHb). 
The insets of FigEIa,b) show sketches of the lateral and 
vertical views of two cylindrical dots with vectors 6, d 
pointing from the smaller dot to the larger one. The 
vector (3 is oriented in the direction 6xd, and for small 
differences Aa=a + — a_ between the dot radii a±, it is 
proportional to Aa. It follows from Eq.QJ that for two 
cylindrical dots of equal sizes (3/ J vanishes because the 
system then has a center of inversion and the ground state 
is symmetric is no asymmetric exchange because the sys- 
tem has cylindrical symmetry (6xd=0). In each of these 
figures there are two regions of behavior: (i) For small 
size differences Aa (the left sides of the peaks) (3/ J in- 
creases roughly proportional to d in FigQf a) and roughly 
proportional to b in FigEJb) until it reaches ~3.5xl0 -3 ; 
in this region the two electrons are distributed almost 
symmetrically between the dots, (ii) For larger values of 
Aa (right sides of the peaks), (3/ J first increases with d 
to a small maximum and then decreases [FigEJa)], and 
similarly with b Fig^b)] . In this second region both elec- 
trons tend to occupy preferentially the larger dot, where 
the Coulomb energy is overcome by the difference be- 
tween single-particle energies of the two dots. 

In order to study shape asymmetry, we consider two 
vertically coupled dots with deviations from cylindrical 
symmetry. In FigOwe give results for two identical el- 
liptical dots that are rotated by tt/2 with respect to one 
another, as shown in the inset, with an offset 6=4 nm. 
The dots of equal sizes have equal energies, which leads to 
an equal distribution of two electrons on them. The ratio 
(3/ J reaches substantial values (~10 -3 ) as a function of 
the separation d and has a relatively weak dependence on 
the angle 7 between the axis of the dots. In this case (3 
arises from the shape asymmetry. The coupling depends 
on the angle 7 between the relative position vector and 
the principal axis (for cylindrical dots, (3 changes direc- 
tion but is constant in magnitude). The orientation of (3 
is given by bxd. In the case where the axies of the dots 
are parallel (a-^ x =a+^ x , a_ 52y =a +52y ), the asymmetric ex- 
change is zero (the system then has an inversion center 
at (6/2, d/2)). In cases when the dots are different in size 
and not cylindrical, (3/ J is even larger. For example, for 




FIG. 2: (Color online) Dependence of modulus of asymmet- 
ric exchange (3 /J on angle 7 between the axes of two ellip- 
tical vertically-coupled dots of equal major and minor axes 
a_ )CC =a+,^=3.5 nm, a_, y =a+, x =6.5 nm, at b—A nm. The in- 
set gives a projection of the dot contours on a plane per- 
pendicular to the growth axis, using the harmonic potential 
parameters to represent the dot sizes. 



a cylindrical dot with a_=5 nm coupled with an elliptical 
dot with a +5 ^=4 nm, a +52y =6.25 nm, the maxima of (3/ J 
are ~3.5xl0 -3 for all 7. 




FIG. 3: (Color online) Dependence of modulus of asym- 
metric exchange (3 /J on angle 7 between the axes of 
two elliptical laterally-coupled dots of equal potential sizes 
A-, x =A+, y =6 nm and A-, y =A+, x =9 nm. The inset shows 
the two QDs that lie in the same plane, seen along the growth 
axis, using the Gaussian potential parameters to define their 
sizes. 

We have calculated the contribution to the asymmetric 
exchange for these structures from the bulk Dresselhaus 
coupling H^ ) =i r y^ ) d x (dy — d^)S a /h (plus cyclic permuta- 
tions of cartesian indices). The coupling constant 7^ is 
47 meVnm 3 for GaAs, and of the order of 100 meVnm 3 
for InAs. In FigQJc) we compare the contribution (3d/ J 
to the asymmetric exchange from the Dresselhaus cou- 
pling with asymmetric exchange (3/ J for a size difference 
Aa=0.22 nm |H El & We see that P/ J is lar £ er for 
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intermediate separations d, a region of particular inter- 
est for implementations for quantum information. As the 
difference in the dot sizes increases, (3/ J becomes larger 
relative to /3d/ J •> and as the size difference decreases, 
(3d/ J becomes larger. 

We have also considered laterally coupled QDs. To ob- 
tain a representation of the barrier between the QDs we 
use inverted-gaussian potentials [22j], and we use again 
the material parameters for In As. The wavefunctions 
are obtained variationally [2?j. In Fig.© we give results 
for two elliptical dots of equal sizes rotated with respect 
to each other by it/ 2. In this case the anisotropic ex- 
change (3/ J arises from the shape asymmetry. From the 
operator rxVi in Eqd the asymmetric exchange f3 has 
a nonzero component only along the growth axis, and 
the dependence of the modulus (3/ J is symmetric with 
respect to 7=7r/4. In Fig.(|3J) once again (3/ J reaches a 



maximum of ^10 -3 . In cases when the Dresselhaus [27) 
and Rashba [2^ couplings have equal coupling constants, 
their contribution is small for 7=7r/4 and then the total 
asymmetric exchange is dominated by (3/ J. 

In summary, we derived an asymmetric contribution 
to the exchange interaction between two electrons in III- 
V semiconductors that arises from the Coulomb inter- 
action and the band mixing and does not require inver- 
sion asymmetry. For asymmetric coupled semiconductor 
QDs, this contribution depends on the geometry and is 
typically 10 -3 of the isotropic exchange J. This interac- 
tion also can play a role in the relaxation and dephasing 
of spin in transport processes in low-dimensional struc- 
tures. 
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